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Abstract 



The parton spectra as predicted by the ARIADNE Monte Carlo generator, for both e^e annihilation 
and deep inelastic scattering, are compared to the QCD MLLA calculations. 
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1 Introduction 

The perturbative QCD approach to describing the inclusive energy spectra, via the modified leading log ap- 
proximation (MLLA) in conjunction with local parton hadron duality (LPHD), has been very successful in both 
e+e~ annihilation and deep inelastic scattering experiments Q. Using LPHD, the non-perturbative eflPects of 
such distributions are reduced to a simple factor of normalisation that relates the hadronic distributions to the 
partonic ones. Perturbative features of these distributions are calculated by MLLA which accounts for both the 
double and single logarithmic effects. The MLLA approach has two free parameters: a running strong coupling, 
governed by a QCD scale A, and an energy cut-off, Qo, below which the parton evolution is truncated. 

The MLLA evolution equation allow the parton spectra for the logarithmic scaled energy spectra, ^, to be 
calculated The variable ^ is defined as ln(£'o/-E) = \n{l/xp), where Eq is the original energy of the jet and 
E is the parton's energy. The cut-off, Qo, bounds the parton energy, E > kx > Qo, where fcr is the transverse 
energy of the decay products in the jet evolution. In order to reconstruct the ^ distributions one has to perform 
the inverse Mellin transformation: 

m,Y,\)^ rt--z?(c.,r,A) (1) 

where the integral runs parallel to the imaginary axis on the right of all singularities in the complex w— plane, 
Y ^ HEa/Qo) and A = ln(Qo/A). 

The Mellin-transformed distributions, D{uj,Y, X), can be expressed j2| in terms of confluent hypergeometric 
functions, $: 

Diu, y. A) = gf^H-A + B + hB + 2; -h^A -B,l-B; h) 
+ (if) H-A, -B- -hMA, B + l-M), 

where 

ti = uj{Y + A) , t2= oj\. 

In addition A and B are defined as: 



(2) 
(3) 



A = mjbuj, B = a/b, (4) 
where Nc is the number of colours, a = llNc/'3+2nf /3N^ , Uf is the number of fiavours and b = llNc/3+2nf /3. 
Equation |^ is then calculated using a numerical integration in the complex plane. 

The current region in the ep Breit frame is analogous to a single hemisphere of e+e^ annihilation. In 
e+e~ — > qq annihilation the two quarks are produced with equal and opposite momenta, ±y^/2, where ^/s is 
the positron-electron centre of mass energy. The fragmentation of these quarks can be compared to that of the 
quark struck from the proton; this quark has an outgoing momentum ~Q/2 in the Breit frame, where is the 
negative square of the four- momentum of the virtual exchanged boson in DIS. In the direction of this struck 
quark the scaled momentum spectra of the particles are expected by MLLA |^, || to have a dependence on 
Q similar to that observed in e"*"e~ annihilation at energy ^/s = Q, with no Bjorken-x dependence. 

The ARIADNE Monte Carlo generator §| is based on the colour dipole model, CDM 0. In the CDM, aU 
gluon emissions constituting the QCD cascade start as radiation from the colour dipole formed between the 
quark and the anti-quark in the case of e~^e~ annihilation or the struck quark and the proton remnant in the 
case of DIS. All subsequent radiation arises from independent colour dipoles formed either from qq pairs or 
softer gluons radiated by the previously produced gluons. In the DIS scenario, the proton remnant is treated as 
an extended object which results in a suppression of radiation generally in the proton direction. In addition 
the struck quark is treated as extended, as the photon only probes it to a distance inversely proportional to the 
transferred momentum. Treating the remnant and the struck quark as extended objects, rather than point like, 
results in a reduction in the available phase space for gluon radiation in DIS. 

The QCD cascade in ARIADNE is governed by a number of parameters in the Monte Carlo models. Two 
of the most important are the QCD scale. A, (PARA(l)) and the parameter that determines the kr cut off for 
the shower (PARA(3)). An additional parameter, PARA(28), also allows the user to bound the lower energy of 
the emitted parton as well. For this study ARIADNE version 4.10 has been used. 



2 Comparisons with MLLA 



Before investigating the evolution of the shower in DIS, the evohition in the simpler case of qq pair production 
in e+e" annihilation was studied. The spectra for both MLLA and ARIADNE were generated with a A = 150 
MeV and a cut-off Qo = 2A = PARA(3) = PARA(28). Below ^ w 1 there are instabilities in the numerical 
integration of equation ^ so all subsequent comparisons are for ^ > 1. Except for an overall normalisation 
discrepancy (a factor of 1.4 greater parton multiplicity in ARIADNE) the ^ spectra are in very good agreement 
as illustrated in Figure |^. This normalisation discrepancy is constant, independent of the ^/s at which the 
events were generated at. There is a slight tendency for the MLLA calculation to fall off quicker at large values 
of £, than the ARIADNE predictions. 
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Table 1: The {x,Q) analysis bins. 

Using LEPTO to generate the electroweak cross section and colour flow configuration for DIS, ARIADNE 
was then used to generate the subsequent QCD cascade. The event was boosted to the Breit frame and those 
partons in the current fragmentation region selected. The DIS events were generated with fixed kinematics that 
are accessible in the HERA regime. The corresponding value of Bjorken-x with Q are shown in Table [|. Using 
the same values of A and Qo that was used for the e+e^ annihilation study, the MLLA prediction is again 
compared to the ARIADNE generated spectra. 

Figure H shows the default version of ARIADNE for DIS compared to the MLLA predictions. As increases 
the discrepancy between ARIADNE and the MLLA calculations becomes more pronounced. The ^ distribution 
of ARIADNE peaks at higher values than the MLLA calculation. In addition, the MLLA calculations are 
narrower than the ARIADNE predictions. Again the parton multiplicity of the two distributions are different. 
Unlike the e+e^ situatuion this normalisation factor seems to exhibit a Q dependence. At low Q the height of 
the peak for ARIADNE compared to MLLA is a factor of 1.1 higher whilst in the highest Q bin it is 1.1 lower. 

In the default ARIADNE, the mechanism for soft suppression of radiation due to the extended source of the 
proton remnant results in a suppression of radiation in the current region of the Breit frame at high Q^. Figure ^ 
shows the high modified version of ARIADNE |]l^ for DIS, where this suppression in the current region is 
removed, compared to the MLLA predictions. As expected, this modification to ARIADNE leads to a much 
better agreement between the MLLA calculations and ARIADNE. The situation with the parton multiplcity is 
similar. The Q— dependence of the ratio of the peak heights is less than the default ARIADNE, with ARIADNE 
being a factor 1.2 — 1.3 higher. 

In both options of the ARIADNE program there are discrepancies evident in the lower (x, Q) bins compared 
to the MLLA predictions. One possible cxplaination of this discrepancy is given in Ref. ||ll|, where it shown that 
high pt emissions in DIS can lead to the situation where the current region of the Breit frame is depopulated. 

3 Conclusions 

The QCD cascade as implemented in the ARIADNE Monte Carlo program is in good agreement with the shape 
of the MLLA prediction in the simple scenario of qq production in e^e~ annihilation. In the more complex 
situation of DIS the agreement is not as good, unless account is taken of the additional suppression introduced 
into the model in the current fragmentation region caused by the suppression of phase space due to the extended 
nature of the proton remnant. 
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Figure 1: The logarithmic scaled momentum distributions, ^, for the ARIADNE Monte Carlo generator (full 
line) and MLLA calculation (dashed line) for e+e" annihilation. The ARIADNE prediction has been scaled 
down by a factor 1.4. 





Figure 2: The logarithmic scaled momentum distributions, ^, for the default ARIADNE Monte Carlo generator 
(full line) and MLLA calculation (dashed line) for DIS. The distributions have been normalised to the peak 
height of the MLLA calculation. 




Figure 3: The logarithmic scaled momentum distributions, ^, for the modified ARIADNE Monte Carlo gen- 
erator (full line) and MLLA calculation (dashed line) for DIS. The distributions have been normalised to the 
peak height of the MLLA calculation. 



